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Controlled manipulation of quantum states is central to studying natural and 
artificial quantum systems. If a quantum system consists of interacting sub- 
units, the nature of the coupling may lead to quantum levels with degenerate 
energy differences. This degeneracy makes frequency-selective quantum oper- 
ations impossible. For the prominent group of transversely coupled two-level 
systems, i.e. qubits, we introduce a method to selectively suppress one transi- 
tion of a degenerate pair while coherently exciting the other, effectively creating 
artificial selection rules. It requires driving two qubits simultaneously with the 
same frequency and specified relative amplitude and phase. We demonstrate our 
method on a pair of superconducting fiux qubits [1]. It can directly be applied to 
the other superconducting qubits p"^? and to any other qubit type that allows 
for individual driving. Our results provide a single-pulse controlled-NOT gate 
for the class of transversely coupled qubits. 



Transverse coupling implies an interaction that is transverse to the eigenstates of the 
uncoupled systems. For this type of coupling the energy splitting of one qubit does not 
depend on the state of the other. This property is appealing because it means that in 
the absence of driving the system essentially behaves as a set of uncoupled qubits. The 
coupling primarily manifests itself when the system is driven and hence can be regarded as 
AC-tunable [7|. A coupling that is not purely transverse leads to a spectroscopic splitting 
of the transitions. Although this splitting enables simple resonant driving for all operations 



91], in practice it requires refocussing schemes to compensate for the continuously evolving 



phases 



10|. The price to pay for the advantage of transverse coupling is obvious; the 



degeneracy prohibits schemes for selective excitation that rely on a frequency splitting. 
Previous experiments used either additional coupling elements [ill], extra modes 12], or 
shifted levels into and out of resonance by DC [13] or strong AC fields 15]. Note 
that level shifting can imply passing through conditions of low coherence 16|], or passing 
resonances with other qubits. In contrast, our method works for simple direct coupling as 
well as for systems with additional coupling elements, such as harmonic oscillators, as long 
as the effective coupling is transverse. It uses only a single pulse of a single frequency and 
does not require (dynamical) shifting of the levels. 
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We consider the class of systems of transversely coupled qubits, described with the Hamil- 
tonian 

H = {A,al + A,al) + Jalal (1) 

where Aj is the single-qubit energy splitting of qubit i, Ai 7^ A2, J is the qubit-qubit 
coupling energy and crj.^^ are t 



actively used quantum systems 



le Pauli spin matrices. This Hamiltonian describes many 
l|-|6|], and often applies for operation at a coherence sweet- 
spot [3I. I14J. Il7l4l9l|. The energy levels of this system are shown schematically in Fig. lb. The 
arrows indicate the transitions of interest; the blue and red arrows describe the transitions 
of qubit 1 and 2, respectively. Both pairs are degenerate in frequency, which is typical for 
transverse coupling. 

Our method aims at the selective excitation of these transitions and is based on simultane- 
ously driving the two qubits with a single resonant frequency, employing different amplitudes 
and phases. The driving is described with the Hamiltonian 

Hdrive= aicos{ut + (pi)al. + a2Cos{ut + (p2)al, (2) 

where u is the driving frequency, and and are the driving amplitude and phase for qubit 
i. The transition strength Tk^i = (/|ifdrive|^), with the driving Hamiltonian transformed to 
an appropriate rotating frame (see Supplementary Information), governs the transition rate 
and depends on both and ipi. Figure Id shows the normalized |T| = \T\/{ai + 02) as a 
function of 01/(01+02), for a fixed phase difference ip2—'Pi = 0. Clearly the two transitions of 
each qubit generally do not have the same strength T, despite their frequency degeneracy. In 
addition, for certain settings individual transitions are completely suppressed: the transition 
is darkened (black dotted lines in Fig. Id). The darkened transitions provide the desired 
conditions where one of the two transitions can be excited individually, even though the 
driving field is resonant with both transitions. The difference in transition strength can be 
understood intuitively. As coupling leads to mixing of the single-qubit eigenstates, qubit 2 
can be excited by driving qubit 1 with a frequency that is resonant with qubit 2. This indirect 
driving can enhance, counteract and even cancel direct driving of qubit 2. The effect differs 
for the two degenerate transitions, because the states involved are different superpositions 
of the single-qubit eigenstates (see Supplementary Information). For J ^ |Ai — A2I, and 
assuming Ai > A2, one readily finds the driving amplitude ratio 

0-2 J /„N 



which yields Tqo^oi = for (^92 — V^i = and Tio^n = for — ^\ = vr. These are the 
transitions of qubit 2. For the transitions of qubit 1, i.e. 00 o 10 and 01 o 11, the amph- 
tude ratio is simply inverted. Expressions for arbitrary J are given in the Supplementary 
Information. 

To experimentally demonstrate this method we employ two coupled flux qubits [l|, each 
consisting of a superconducting loop interrupted by four Josephson tunnel junctions. When 
biased with a magnetic flux close to half a flux quantum $0; the two states of each qubit are 
clockwise and anti-clockwise persistent-current states. These currents Jp produce opposite 
magnetic fields, which provides the coupling for the two qubits. Two independent AC- 
operated SQUID magnetometers are used to simultaneously read out the states of the qubits 



20I, 



21| . These are switching-type detectors, where the switching probability Pg^ is a measure 



for the magnetic field. At a bias of ^ the system is described by the Hamiltonian of Eq. 1. 
Here the eigenstates of each qubit are symmetric and anti-symmetric superpositions of the 
two persistent-current states, with level separation A. The device is depicted in Fig. la. 
The qubits are characterised by the persistent currents Jp^i = 355 nA and /p,2 = 460 nA and 
the energy splittings Ai//i = 7.88 GHz and A2//i = 4.89 GHz. The qubit-qubit coupling 
strength is 2J//i = 410 MHz. 

For our fabricated quantum objects the spatial locations are well-defined, and the in- 
dividual control of amplitude and phase for each qubit according to Eq. 2 can be easily 
achieved using local magnetic fields. We employ two on-chip antennas, indicated as A\ and 
A2 in Fig. la, both coupling to both qubits, with a stronger coupling to the closer one. 
Driving the two qubits from both antennas is described with 

ifdrive = M COs{ut + 0i)(miiCT^ + mi2cr^) 

+A2 COs(wt + 02) ("^2lO-x. + ^22(^1) , (4) 

where Aj and (pj are the driving amplitude and phase for antenna j and rriji is the coupling 
of antenna j to qubit i. Note that any combination of and ipi in Eq. 2 can be achieved 
with the proper choice of Aj and (pj. For this device mu/mu = 0.32, m2i/m22 = 0.33 and 
mil = "^22- 

For the experimental demonstration we choose to focus on the degenerate transitions of 
qubit 2. We first show that, if the qubits are driven from a single antenna, the two degenerate 
transitions exhibit a different Rabi frequency. We apply two pulses on antenna 1: the first 



pulse is resonant with qubit 1, the second pulse is resonant with qubit 2. Figure Ic shows 
a schematic of the pulse-sequence; note that here A2 — 0. The experiment is repeated for 
varying durations Ti and T2 of pulses 1 and 2. The switching probability Psw,i of detector 
1 is depicted in Fig. 2a, showing a few Rabi oscillation periods as a function of the pulse 
duration Ti. Varying T2 does not lead to oscillations of qubit 1, as pulse 2 is non-resonant, 
and only relaxation is observed. The oscillations of qubit 2, induced by the second pulse, 
are visible in Psw,2 (Fig. 2b). Here we distinguish two oscillation frequencies. Along the 
white solid line, where qubit 1 is prepared in the excited state, qubit 2 oscillates with a Rabi 
frequency of 85 MHz. For qubit 1 prepared in the ground state, along the white dashed line, 
the Rabi frequency / = 31 MHz of qubit 2 is lower. For qubit 1 in a superposition of the 
ground and excited states, qubit 2 shows a beating pattern of both oscillations. 

A more detailed analysis allows us to unravel the two frequencies of Fig. 2b and determine 
which levels are participating in each of the oscillations. We extract the level occupations 
Qoo, Qoi, Qio and Qu from the individual switching probabilities (see Supplementary Infor- 
mation). The result is shown in Figures 2c- f. After an odd number of 7r-rotations of qubit 
1, there are only oscillations between states 10 and 11, not for states 00 and 01. After an 
even number of vr-rotations of qubit 1 the situation is reversed; now the states 00 and 01 
oscillate. The two oscillation frequencies are clearly linked to the two different transitions. 

To demonstrate the tunability of the transition strengths we drive both antennas simul- 
taneously, using the same frequency and controlling independently the amplitudes ^1,^2 
and phases 0i,02- In this 2-pulse experiment (Fig. Ic), the first pulse prepares qubit 1 
with a 7r/2- rotation and the duration T2 of the second pulse is varied. Since qubit 1 is in a 
superposition state, both Rabi frequencies are present in the dynamics of qubit 2. In Figures 
3a-c we show the Fourier transform for the measured oscillations. Each graph is measured 
with a different amplitude ratio A1/A2, with fixed phase (f>i — and varying 02- Figure 3a, 
with A1/A2 — 1.3, shows a typical result for an arbitrary amplitude-ratio; the Rabi oscil- 
lation frequencies of both transitions clearly depend on 02 — 0i, but nowhere a transition 
is darkened. Note the occurrence of equal Rabi frequencies for two phase conditions, as 
denoted by Y. For A1/A2 = 2.5 in Figure 3b we observe that for 02 — 0i ~ tt (indicated 
by Xq) the transition 10 -H- 11 is fully darkened, while the 00 ^ 01 transition shows a 
non-zero oscillation frequency. In Fig. 3c with A1/A2 — 6.3 the situation is reversed, with 
the 00 •<-> 01 transition being suppressed (denoted by Xi). This clearly demonstrates our 
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method, as we selectively excite one of two transitions, despite their frequency degeneracy. 
Calculations of T are in good agreement with the experimental results, provided we allow 
for different transmissions of amplitudes and phases of the antennas to the qubits, which we 
attribute to the influence of the detector circuits. 

To further investigate the special cases of equal Rabi frequencies (Y), and darkened tran- 
sitions (Xo, Xi), we again vary the durations ti and T2, using both antennas for the second 
pulse. The results should be compared with Fig. 2b. Figure 3d shows Psw,2 for driving 
conditions denoted by Y (left arrow): the oscillation frequency of qubit 2 does not depend 
on the state of qubit 1. For the conditions marked by Xq, we only observe oscillations of 
qubit 2 when qubit 1 is in the ground state, as shown in Fig. 3c. Similarly for the conditions 
marked by Xi, now we only observe oscillations of qubit 2 if qubit 1 is in the excited state 
(Fig. 3f). 

The demonstrated capability to selectively manipulate transition strengths in frequency- 
degenerate transitions has important applications. A 7r-pulse using condition Xi or Xq 
provides a l-controUed and O-controUed NOT gate, respectively. This enables certain sys- 
tems, including the flux qubit used here, to be fully operated at the coherence-optimal 
point, without level shifting by either DC or strong AC signals. Note that the use of addi- 
tional coupling elements is neither required nor prohibited. If additional coupling elements 
are used, our method can replace more complicated schemes. For conditions similar to Y , 
taking care of the individual rotation angles, also single-qubit gates can be implemented. 
The controlled- NOT and single-qubit gates together form a universal set, implying that our 
method fulfills all requirements for constructing any single- or two-qubit gate. The method 
also scales to three or more qubits, provided that for a certain target pair the system can 
be reduced to Eq. 1. 

In conclusion, we have introduced and experimentally demonstrated a method to con- 
trol transition strengths by applying a non-uniform driving field. Darkened transitions are 
created and employed for the selective excitation of degenerate transitions. As this method 
improves the simplicity and coherence conditions for operations in a variety of quantum sys- 
tems, the prospect of performing large-scale quantum algorithms is enhanced significantly. 
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FIG. 1: Coupled qubit system and transitions, a, Optical micrograph of the sample, showing 
two flux qubits colored in blue and red. The inset shows part of each qubit loop, both containing 
four Josephson tunnel junctions. Overlapping the qubit loops, in light-grey, are the SQUID-based 
qubit-state detectors. In the top right and bottom left are the two antennas from which the qubits 
are driven, b. Energy level diagram of the coupled qubit system. Arrows of the same color indicate 
transitions of the same qubit and are degenerate in frequency, c. Pulse sequence used for the 
coherent excitation of the qubits. The first pulse is resonant with qubit 1. The second pulse, applied 
from both antennas simultaneously with independent amplitudes and phases, is resonant with 
qubit 2. After the second pulse the state of both qubits is read out. d. The normalized transition 
strengths of the four transitions in b as a function of the net driving amplitudes ai/(ai + 02) for 
(p2 — ^1 = 0. For ip2 — ifi = the dashed and solid lines are interchanged. The black dotted lines 
indicate the locations of the darkened transitions. 
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a b 




pulse 1 duration t, (ns) 

FIG. 2: Driving from a single antenna. Measurement of the state of the qubits, represented 
by switching probabilities Psw,i and Psw,2, after applying a pulse of duration ti resonant with qubit 
1, followed by a pulse of duration T2 resonant with qubit 2. a, Psw,i, showing coherent oscillations 
of qubit 1 induced by pulse 1. The white solid and dashed lines indicate a tt- and 27r-rotation 
respectively. For pulse 2, qubit 1 only shows relaxation, b, Psw,2, showing coherent oscillations 
induced by pulse 2. After an odd number of 7r-rotations on qubit 1, the oscillation frequency is 
higher than after an even number of vr-rotations. For superposition states of qubit 1, a beating 
pattern of the two oscillations is observed, c-f, Level occupations Q of the four different levels. 
Note that a value of 0.2 has been added to Qu to improve visibility. 
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FIG. 3: Transition strength tuning and darkened transitions, a-c, Rabi frequency de- 
pendence on (j)2 — (pi for three different amplitude-ratios. The color scale represents the Fourier 
component of Psw,2iT2)- Qubit 1 is prepared with a 7r/2-rotation. Markers Xq and Xi indicate the 
conditions for a darkened transition on 00 o 01 and 10 -H- 11 respectively, d-f Psw,2 versus the 
durations ri and T2. The white solid and dashed lines indicate a vr- and 27r-rotation of qubit 1, 
respectively. The driving conditions are as marked by Y left arrow (d), Xq (e) and Xi (f). 
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Supplementary Information 



CONDITIONS FOR DARKENED TRANSITION 



Let us take the Hamiltonian: 

if = 



(1) 



with Ai 7^ A2. Without loss of generahty we take Ai > A2. Additionally we have the 
driving term: 



drive 



ai cos{ujt + ^i)crl + 02 cos(a;t + ^2)crl 



H. 



drive 



H. 



drive 



where 



drive 



-iu)t 



a2 



2' " ■ 2 

The Hamiltonian (Eq. [1]) can be diagonalized without any approximations: 



|0) 
|1) 
|2) 
|3) 



cos^ilOO) -sin^illl) 
cos ^2 101) - sin ^2 1 10) 
cos ^2 1 10) + sin ^2 1 01) 
^i|ll) + sin^i|00), 



cos 



(2) 
(3) 

(4) 

(5) 
(6) 
(7) 
(8) 



where 



(9) 
(10) 



2J 

tan 2^1 = — 

A1 + A2 

2J 

tan 26*2 = • 

Ai - A2 

To calculate the transition strength T^^^i = (/|ifdrive|^), we ignore the counter- rotating fields, 
which results in i/drive = -f^drive for /c < / and i^drive = -^drive ioT k > I. We find that 

(l|^drive|0) = -ye*'^' {cos^isin02 + sin^iCos^2} + ye^'^^ {cos^i cos02 + sin^i sin02} 

(11) 



2 

= -ye^'^^ sin(ei + 



O2) + y e^^^ cos(^2 



(3|^drive|2) = +ye*'^^sin( 



'1 + ^2) + y e^'^^ cos(^^2-^^i), 



(12) 
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and 

(2|^drive|0) = ^e^^i {cos 9^ cos 62 - sin 9^ sin 62} + ^e^'^^ {cos Oi sin ^2 - sin Oi cos ^2} 

= y e^"^! cos(^i + 82) + y e^'^^ sin(^2 - ^1) (13) 

(3|^drive|l) = ye^^i cos(^i + 62) - ye*^^ sin(^2 - ^1). (14) 

We therefore find that the required ratio of driving amplitudes to drive only one of the two 
transitions of qubit 2 is given by 

a2 _ sin(6'i + 62) , . 

ai ~ cos(^2 - ^1) ' ^ ' 

with (/72 — = for suppressing the 00 ■<->■ 01 transition, and ^2 — ^\ = for suppressing 
the 10 ^ 11 transition. For driving the transitions of qubit 1, the ratio is 

a2 _ cos(6'i + 62) , . 

ai ~ sin(^2 - ^1) ' ^ ' 

with (/?2 — = TT for suppressing the 00 •<->■ 10 transition, and (/72 — = for suppressing 
the 01 11 transition. 

The resonance frequencies are given by: 

E,-Eo^E-,-E2 = ^ (^4 J2 + (Ai + A2)2 - ^4J2 + (Ai - A2)2) (17) 
E2-Eo^E-i-E, ^ ^ (y4J^T(ArTA^+ y4J^T(Ar^A^) . (18) 

If we now make the assumption that J <C |Ai — A2I, then to lowest order approximation 
^1 = and 02 — J/\A.i — A2\. The eigenvectors reduce to 

|0) = |00) (19) 
ID = 101) -^^110) (20) 

|2) = |io) + a73a;I°^^ (21) 

|3) = 111), (22) 

and El — Eo — E3 — E2 — A2, E2 — Eo — E3 — El — Ai. FoUowing the same procedure as 
above we find for the darkened transitions 

- = (23) 
a2 A1-A2 ^ ^ 
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and 

- = -X-^ (24) 

for suppressing the transitions of qubit 1 and 2, respectively. The phase conditions are the 
same as above. 



II. EXTRACTING LEVEL OCCUPATIONS FROM DETECTOR 
SWITCHING-PROBABILITIES 



The detectors arc click-type detectors; they cither switch or do not switch away from the 
initial state, depending on the sensed magnetic field. Their fidelities are lower than 1, and 
the base levels for the switching rates are free to choose. Here we present a procedure to 
extract the level occupations Qi — \ {i\'4>)\'^, with i — 00, 01, 10, 11 and \'4>) the quantum state 
of the system. 

Each individual measurement event can have one of four possible outcomes: neither 
detectors switches, only detector 1 switches, only detector 2 switches, or both switch. We 
determine the respective probabilities PqO) -Poii -Pio and Pu by repeating the measurement 
many times. The individual switching probabilities of the detectors are calculated from the 
combined probabilities: Psw,i — Pw + Pu and Psw,2 — Poi + Pu- 

The quantities Qn and P„ are related by 



P = MQ 



(25) 
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V 



^00 -foi 


-Pio Pu ] 
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^2(1 


^1(1 
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^1(1 


-E2) 
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^1(1 
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£'1^2 





-E2) 

-E,) 
-E2) 
E1E2 



\ 



(26) 



The parameters Gi and Ei arc the switching probabilities for detector i when the corre- 
sponding qubit is in the ground state (G) or the excited state {E) respectively. Assuming 
no crosstalk between the detectors, these four parameters fully characterize the measure- 
ment. The calibration values are determined from independent measurements. Measuring 
the switching probability without applying any driving pulse provides Gi — 0.507 and 
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(j2 = 0.487. In another two experiments either qubit 1 or qubit 2 was resonantly excited 
with a pulse duration much longer than the coherence times of the qubits, so that the fi- 
nal state is a 50/50 incoherent mixture of the ground and excited state. The full excited 
state must have twice this signal. This procedure is more reliable than doing a coherent 
vr rotation, since that approach is susceptible to gate errors. We measure Ei = 0.760 and 
E2 = 0.744. 

The matrix M is invertible for Gi ^ Ei, and the equation 

Q = M-^P (27) 

is used to extract the desired level occupations. Supplementary Figure 1 shows the data for 
the four measurement outcomes of the experiment described in Fig. 2. 




pulse 1 duration (ns) 

FIG. SI: a-d Measured combined probabilities P in color versus the durations ri and T2 of pulse 
1 and pulse 2. 
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